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Aggregation phenomena are ubiquitous in nature, encompassing out-of-equilibrium 
processes of fractal pattern formation, important in many areas of science and tech¬ 
nology mm- Despite their simplicity, foundational models such as diffusion-limited 
aggregation jSj (DLA) or ballistic aggregation [6j (BA), have contributed to reveal the 
most basic mechanisms that give origin to fractal structures. Hitherto, it has been 
commonly accepted that, in the absence of long-range particle-cluster interactions, the 
trajectories of aggregating particles, carrying the entropic information of the grow¬ 
ing medium, are the main elements of the aggregation dynamics that determine the 
fractality and morphology of the aggregates [71-fTT]. However, when interactions are 
not negligible, fractality is enhanced or emerges from the screening effects generated 
by the aggregated particles, a fact that has led to believe that the main contribution 
to fractality and morphology is of an energetic character only, turning the entropic 
one of no special significance, to be considered just as an intrinsic stochastic ele¬ 
ment [121420] . Here we show that, even when long-range attractive interactions are 
considered, not only screening effects but also, in a very significant manner, parti¬ 
cle trajectories themselves are the two fundamental ingredients that give rise to the 
fractality in aggregates. We found that, while the local morphology of the aggre¬ 
gates is determined by the interactions, their global aspect will exclusively depend 
on the particle trajectories. Thus, by considering an effective aggregation range. A, 
we obtain a wide and versatile generalization of the DLA and BA models. Further¬ 
more, for the first time, we show how to generate a vast richness of natural-looking 
branching clusters with any prescribed fractal dimension D G [l,Do], very precisely 
controlled, being Dq the fractal dimension of clusters resulting from the underlying 
non-interactive aggregation-model, DLA or BA, used. 
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Since the introduction of the DLA model, a plethora of studies and models have been 
developed, trying to understand the most basic mechanisms that give rise to self-affine 
structures in aggregation phenomena M- In particle-cluster aggregation, the paradig¬ 
matic DLA and BA models have shown that, if long-range interactions are negligible, the 
characteristics of the growing medium, such as its temperature and viscosity, encoded in 
the aggregating particle trajectories, will dehne the morphology and fractality of the grow¬ 
ing cluster [Hi. This fact is remarkably exhibited when, by changing the particles mean 
squared displacement (MSD), a transition from compact BA clusters to branching DLA ones 
can be achieved [SHn]. On the other hand, when interactions can no longer be disregarded, 
aggregation dynamics could become quite complex. Nonetheless, experimental and compu¬ 
tational models have shown that the effects of interactions on the morphology and fractality 
of clusters are mainly two. Short-range repulsive interactions decrease the fractality, as they 
allow particles to reach a minimum in the energetic landscape, thus producing compact 
clusters [T^ [TH] . Conversely, long-range attractive interactions increase fractality, gener¬ 
ating highly ramihed clusters, due to screening effects coming from the interaction-range 
of the aggregated particles [la EHm EiHS]. However, one must carefully consider that 
the main mechanism responsible for the cluster’s morphological and fractal features cannot 
trivially be attributed to a single process of an energetic or entropic nature, when the un¬ 
derlying aggregation dynamics comes up as a result of the interplay between the stochastic 
trajectories of the particles and the quasi-deterministic particle-cluster interactions [m ITB] . 
While several technically and conceptually simple approaches have regarded the energetic 
character of the dynamics as the most important factor to the clusters morphology and 
fractality [IMllIIMQ], in this Letter we show that the entropic contribution cannot be 
trivially considered as an intrinsic stochastic element of the system, but as an important 
aspect that not only contributes to the fractality of the clusters, but also as a remarkable 
source of diversity in fractal pattern formation. 

To this end, we introduce a two-dimensional particle-cluster model that, trough the incor¬ 
poration of an effective interaction or aggregation range A, deeply generalizes and enhances 
the capabilities of the standard DLA and BA models. By exploring the aggregation dynam¬ 
ics under constant A and A = A(A^), where N is the number of aggregated particles in the 
cluster, this simple but non-trivial stochastic scheme for aggregation allows us to separate 
and characterize the subtle contributions of energetic and entropic character of the dynamics 
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Figure 1. DLA-based multiscaling aggregates, a Aggregates containing N = 150 x 10^ 
particles each, for A = 1,10,100 and 1000 units, visualized at 5%, 10%, 30% and 100% of their total 
size. The blue squares display the multiscaling evolution of the structure, b Radius of gyration, 
Rg, and c fractal dimension, D, versus the number of aggregated-particles, N, in log-log and lin¬ 
ing plots, respectively. Notice that, when A —)> oo, the structure of the aggregates tends to MF 
(D = 1). Each curve was computed as an average over an ensemble of aggregates, d Evolution 
of the growing front for the first two stages of growth (see text), e Typical structure of a ME 
aggregate. 


to the fractality and morphology of the clusters, as well as it reveals features not previously 
seen before [El EH]. The aggregation dynamics of this model is quite simple, assuming that 

(i) each aggregated particle in the cluster is assigned with an interaction radius A, and that 

(ii) as soon as the path of a wandering particle intersects the cluster’s interaction boundary, 
it moves radially and ballistically to the closest particle in the cluster, getting irreversibly 
and rigidly attached. Notice that the interaction boundary is defined by the overlap of the 
individual interaction regions of the aggregated particles (see Fig. Id and 2d). Additionally, 
distance quantities are measured in particle-diameter units, which are fixed to one. In or¬ 
der to characterize the fractal properties of our aggregates, we used the radius of gyration, 
Rg = kN^, where fc is a constant, = 1/D, D being fractal dimension. Typically, D is 
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Figure 2. BA-based multiscaling aggregates, a Aggregates containing N = 300 x 10^ particles 
each, for A = 1,10,100 and 1000 units, visualized at 5%, 10%, 30% and 100% of their total size. 
The blue squares display the multiscaling evolution of the structure, b Radius of gyration, Rg, and 
c fractal dimension, D, versus the number of particles, N, in log-log and lin-log plots, respectively. 
Notice that, when A — )• oo, the structure of the aggregates tends to MF {D = 1). Each curve was 
computed as an average over an ensemble of aggregates, d Evolution of the growing front for the 
first two stages of growth (see text), e Typical structure of a MF aggregate. 

estimated from the slope of plots of Rg vs N in logarithmic scale, or equivalently (as in this 
work), through the numerical derivative of logi?^ with respect to N (see Methods section). 

Let us hrst discuss our results for the case of interactions with a constant range. For 
A = 1, or direct-contact interaction, we have the usual DLA or BA models (see Figs. 1 and 
2, respectively); we will refer to the known fractal dimension of the DLA or BA aggregates 
as Dq = 1.71 or 2, respectively. Further on, when A > 1, two main features will emerge due 
to the interplay of interactions and particle trajectories: a multiscaling branching growth 
and a crossover in fractality, from D = 1 (as A —?• oo) to D = Do (when N —)■ oo). The 
former could be qualitatively characterized by three growth stages. In the hrst stage, when 
the radial size of the cluster is small compared to A, the growth is limited by the interactions. 
The aggregated particles are close enough so that their individual interaction regions are 
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highly overlapped, forming an almost circular envelope around the cluster, and making the 
last aggregated-particles the most probable aggregation points for particles incoming to that 
region. More likely, three arms grow, clearly seen as A —?• oo. This stage is characterized by 
D ^ 1. In the second stage, the envelope starts to develop small deviations from its initially 
circular form, with typically three main elongations or growth instabilities, associated with 
the main branches. When the distance between the tips of two adjacent branches becomes 
of the order of A, a bifurcation process starts, generating multiscaling growth. After that, 
when the interactive envelope develops a branched structure itself (see Fig. Id and 2d), 
particles are able to penetrate into the inner regions of the aggregate and a transition in 
growth dynamics takes place, from interaction-limited to trajectory-limited. In the third 
stage, when the distance among the tips of the main branches becomes much larger than 
A, growth is limited by the mean squared displacement of the wandering particles. This is 
clearly appreciated for A = 100 and 1000 in Figs. 1 and 2. For A = 10, even though the 
same crossover is present, it is masked by the small-cluster-size regime. Notice that in both 
DLA and BA cases, for a larger interaction range, a higher number of aggregated-particles 
is required to recover the asymptotic behavior. Yet, the asymptotic value D = Dq and 
the cluster’s global structure characteristic of DLA or BA will be recovered as N —)■ oo. 
However, as A —)■ oo, the structural and fractal features will tend towards a mean-held (MF) 
behavior [18], characterized by D = 1 and three well dehned branches despite the particle 
trajectories. 

The previous results have two important consequences. Firstly, even though interactions 
leave a strong print in the local cluster’s structure and fractality, the stochastic nature of 
the particle trajectories will ultimately determine their global characteristics. Nonetheless, 
clusters cannot be trivially described by a single fractal dimension as it was previously 
thought naES], since their multiscaling behavior is able to span many orders of magnitude 
in the space occupied by the clusters. Secondly, when long-range attractive interactions are 
introduced in the growth dynamics, the only way to obtain monofractal clusters, with a well 
dehned fractal dimension (smaller than Dq), i.e., with serf-affinity, is to maintain a proper 
balance of energetic and entropic contributions to the fractality. In fact, taking into account 
that the size of the clusters in space is proportional to the radius of gyration Rg oc 
this balance can be achieved by scaling the interaction range with the number of particles 
in the cluster, through the generalized X{N) = XqN^. Here, Aq is hxed to one, while e is the 
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Figure 3. Tunable aggregates. In a and b, log-log plots for Rg and lin-log for D vs N, for 
aggregates grown with specific values of e in the interval [0.01,1], for BA (blue) and DLA (red) 
up to = 10^ particles. Notice how the multiscaling behavior gives way to a single well-defined 
fractal dimension D = D{e). At the bottom, one can appreciate the transition in these monofractal- 
aggregates morphology with respect to e . In c, plots of e vs H for aggregates based on DLA and 
BA. These numerically obtained curves can be used to grow clusters with any prescribed D. d BA- 
and DLA-based clusters are shown with the same fractal dimension, D = 1.51 and 1.31, grown 
with a very high precision around de desired value. 

scaling parameter that takes values in [0,1]. As shown in Fig. 3, this approach precisely and 
uniquely defines D for each given value of e. In addition, using the underlying DLA and 
BA models, within this scheme (scaling A with N)^ allows us to generate a rich variety of 
appealing and novel tunable aggregates with D in [l,iAo], not previously seen [2S1I2Z!- Yet, 
even in this case, the overall morphology of the cluster is still defined by the MSD of the 
wandering particles. 

It is worth to mention that, in contrast with screened-growth or sequential-algorithm 
models [261 Elj, the fractality displayed by the clusters generated with our model is an 
emergent property of the aggregation dynamics. Additionally, our model remains simple in 
contrast to Lagrangian models (i.e., based on a molecular dynamics approach), where the 
cluster fractality comes from the forces (stochastic and deterministic) experienced by the 
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Figure 4. Fractality diagram. The diversity of aggregates that can be obtained with the 
introduction of \{N) is astounding, (i) With e = 0 and Aq = 1, one has the well known transition 
from BA to DLA by changing the trajectories’ fractal dimension, d^, from to 1 to 2, respectively, 
(ii) With e = 0 and a constant Aq > 1 multiscaling aggregates are obtained, while MF behavior is 
obtained in the limit A —?■ oo for a finite N. Otherwise, with a large enough N, we get back to the 
usual BA or DLA behavior. Finally, {in) with e > 0 and Aq = 1, aggregates with a tuned D from 
BA or DLA to MF, can be obtained by adequately scaling A with N. 

particles Moreover, our results allow us to understand, in a very clear manner, the 

effects coming from the wandering-particles’ MSD and from the interactions themselves, on 
the clusters’ morphology and fractality. 

As hnal remarks. Fig. 4 shows the whole family of fractal aggregates that can be gener¬ 
ated under this simple scheme. When the interaction range is kept constant, all multiscal¬ 
ing aggregates will belong to the same universality class corresponding to the underlying 
non-interactive aggregation model used to generate them. On the other hand, when the 
interaction range is properly scaled with the cluster’s size, N, the full range of tunable 
monofractal aggregates can be generated with D G [1,2]. Also notice that our model can 
easily be extended to dimensions higher than 2. Therefore, these very important features of 
our model can be exploited beyond aggregation phenomena. We anticipate our hndings will 
provide important insights into the study of fractal growth phenomena, from networks [28] 
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and branching morphogenesis [29], to bio-inspired materials engineering [30|, among others. 


METHODS 

Attractive interactions. In order to introdnce attractive interactions in onr model, 
an interaction circle of radius A is assigned to each aggregated-particle, centered around 
the particle. This effective interaction distance. A, is measured in particle-diameter units, 
acquiring values A > 1. Additionally, as explained before. It can remain constant all along 
the dynamics or scale with N, depending on the desired application, i.e., to produce a multi¬ 
scaling aggregate or an aggregate with tuned fractal dimensions. Further on, for aggregates 
based on BA, we follow the standard procedure in which particles are launched at random 
from the circumference of a circle of radius L = r^ax + < 5 , with equal probability in position 
and direction of motion m Here, r^ax is the distance of the farthest particle in the cluster 
with respect to the seed particle placed at the origin. In our simulations we used 6 = 1000 
particle diameters to avoid undesirable screening effects. In the case of aggregates based on 
DLA, particles were launched from a circle of radius L = r^ax -l- A -|- 5, with 6 = 100. The 
mean free path for the motion of the particles is then set to one particle diameter, Aq = 1. 
We also used a scheme that modihes (increments) the mean free path as the particles wander 
at a distances greater than L or in between branches, and set a killing radius at Lj. = 2L, 
in order to speed up the aggregation process. In both scenarios, when a particle crosses for 
the hrst time the aggregation zone of any particle, characterized by A, it moves radially and 
ballistically the closest particle along its path in an irreversible (no disaggregation) and rigid 
(no re-accommodation) form, as schematized step-by-step in Extended Data Figure 1. 

Computing the fractal dimension. The fractal dimension, D, is estimated from the 
radius of gyration, Rg, by means of a linear £t of the function Rg = kN^ to the numerical 
data in a log-log plot, where fc is a constant, N is the number of aggregated particles and 
(3 = 1/D. In practice, it is assumed that (3 is constant as long as the number of particles in the 
cluster is large. Because the multiscaling models do not develop a constant fractal dimension 
at their early stages of growth, the simplest way to quantitative and qualitative measure the 
behavior of D is through the derivative of Rg in the logarithmic scale. We did so by means of 
standard two and three point numerical differentiation methods: f'{x) = [/(x-f-h) —/(x)]/h, 
at the ends of the differentiation intervals and f'{x) = [f{x + h)—f[x — h)]/2/i, in between. 


Here f{x) = log Rg{N) and R{g) is computed as a discrete quantity therefore, h is set as 
the distance between the points, x = Xj and x + h = Xj+i. In all cases, Rg is computed as 
an average over a large ensemble of aggregates. Specihcally, the results for the multiscaling 
aggregates (i.e., with A kept constant) shown in Figures la, lb, 2a, and 2b, were obtained 
over 64 and 15 aggregates containing 1.5 x 10^ and 3 x 10^ particles for those based on 
DLA and BA, respectively. In this case, Rg was calculated every 10 particles. In Figures 
Ic and 2c, 192 and 128 aggregates containing 5 x 10*^ and 10^ particles were used to obtain 
de averages, respectively, while Rg was calculated every 7 particles in order to capture the 
features of D at small scales in N. The results for tunable aggregates based on DLA and BA 
(with A = A(A^)), shown in Figure 3, were obtained over 128 and 48 aggregates, respectively, 
containing 10^ particles, and Rg was computed every 10 particles. We must point out that 
the fluctuations observed in D, depicted by the grey curves in Figures Ic, 2c, 3a and 3b, are 
due to the numerical and local aspect of the derivative’s estimation and the stochasticity of 
the model. These fluctuations decrease as N increases, when the aggregates tend towards 
a more well dehned structure. Thus, with the purpose to improve the visualization of the 
observed tendencies, the curves in color were computed by means of a running-average over 
N. 

Tunable aggregates. Aggregates with a prescribed (tuned) fractal dimension, either 
based on BA or DLA, have a well dehned D for each given value of e. Therefore, we computed 
D for e G {0.01, 0.05, 0.1,..., 0.95,1.0} in order to obtain the functional dependence of e on 
D. Then, for a desired D*, we estimated the corresponding value of e through a linear 
approximation using the two closest points in D to the desired D*. In this case, a linear 
approximation is a valid method since the difference among consecutive points in e{D) is 
small and the curve is well behaved as can be appreciated in Figure 4c. 
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EXTENDED DATA 



Extended Data Figure 1. Implementation of the interaction. A step-by-step diagram is 
provided regarding the implementation of the attractive interactions in our model, a Each particle 
in the cluster is assigned an effective radius of aggregation A, starting with the seed particle, b 
Upon aggregation, the independent interaction regions of each particle, defined by A, overlap among 
them, c A particle far from this region does not interact with the cluster until d its trajectory 
is such that, for its next step, it intersects for the first time the interaction boundary of any 
aggregated-particle. This is determined when its perpendicular distance to the particles in the 
cluster is less than A. Then, the position of the aggregated-particles are projected to determine 
the closest one along the direction of motion. Finally, e the position of first crossing is computed 
and the position of the new particle in the cluster is determined. 
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